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Spin-orbit-coupled two-dimensional electron gases (2DEGs) are a textbook example of heli- 
cal Fermi liquids, i.e. quantum liquids in which spin (or pseudospin) and momentum degrees-of- 
freedom at the Fermi surface have a well-defined correlation. Here we study the long- wavelength plas- 
mon dispersion and the Drude weight of archetypical spin-orbit-coupled 2DEGs. We first show that 
these measurable quantities are sensitive to electron-electron interactions due to broken Galileian 
invariance and then discuss in detail why the popular random phase approximation is not capa- 
ble of describing the collective dynamics of these systems even at very long wavelengths. This 
work is focussed on presenting approximate microscopic calculations of these quantities based on 
the minimal theoretical scheme that captures the basic physics correctly, i.e. the time-dependent 
Hartree-Fock approximation. We find that interactions enhance the "plasmon mass" and suppress 
the Drude weight. Our findings can be tested by inelastic light scattering, electron energy loss, and 
far-infrared optical-absorption measurements. 



PACS numbers: 71.45. Gm, 71.10.-w, 71.70.Ej 
I. INTRODUCTION 

In recent years we have witnessed a tremendous explo- 
sion of interest in a large variety of novel two-dimensional 
(2D) quantum many-body systems. Prime examples 
are: (i) strongly spin-orbit-coupled 2D electron and hole 
gases, which are promising candidates for semiconductor 
spintronica^ (ii) graphene 2 ^ (a monolayer of carbon atoms 
arranged in a 2D honeycomb lattice) , which has attracted 
a great deal of interest because of the massless-Dirac- 
fermion character of its carriers and because it may pave 
the way for carbon-based electronics^ (iii) 2D electron 
gases in HgTe/Hg(Cd)Te quantum wells where massless 
Dirac fermions arepredicted to arise at a critical quan- 
tum well thickness 4 ^ and, more recently, (iv) metallic 
surface states of 3D topological insulators^"^. 

These systems share a unique common factor: their 
orbital degrees-of-freedom are intimately coupled to the 
electron spin (or sublattice pseudospin, in the case of 
graphene) degree-of-freedom. This coupling, being of rel- 
ativistic original, naturally breaks Galileian invariance 
and is thus the basic reason for a quite sensitive depen- 
dence of several observables to electron-electron inter- 
actions, even at very long wavelengths (see for example 
Refs.H3 Hl8|) . Furthermore, these systems exhibit coupled 
spin-charge collective dynamics, which is just beginning 
to be investigated in the contemporary literature^. 

In this article we focus our attention on an archetypical 
2D electron gas model Hamiltonian with spin-orbit cou- 
pling (SOC). For the sake of simplicity we choose an ele- 
mentary form of SOC which is linear in momentum and 
has the canonical Rashba or Dresselhaus functional form. 
Since collective dynamics in quantum many-body sys- 



tems is controlled by isolated poles in dynamical linear- 
response functions^, we carry out a microscopic study 
of the density-density response function in the dynami- 
cal limit taking into account exactly SOC and treating 
electron-electron interactions beyond the random phase 
approximation (RPA). The RPA, which is commonly 
used to describe electron liquids, is indeed not capable to 
capture the subtle renormalization of the plasmon mode 
that occurs in non-Galileian-invariant quantum liquids. 
The study of many-body effects when both Rashba and 
Dresselhaus SOC terms are present in the Hamiltonian 
is beyond the scope of the present article: rotational in- 
variance of the Fermi contours is indeed spoiled by the si- 
multaneous presence of both effects and this complicates 
(and partly obscures) the basic interplay between SOC 
and many-body effects we want to highlight. Electron- 
electron interactions in 2D electron and hole gases in the 
pre sence of SOC have attracted a certain deal of atten- 
tio n 1 13 * 14 ^^. Below we will make contact with the pre- 
existing literature whenever possible. 

Our manuscript is organized as follows. In Sect. [TT| we 
present the model we have studied, we introduce the basic 
definitions, and outline the equation-of-motion approach 
we have used to relate the density-density response func- 
tion with the longitudinal current-current response func- 
tion. The latter is then evaluated microscopically within 
the time-dependent Hartree-Fock approximation in the 
long- wavelength limit in Sect. |III| While the main focus 
of this paper is on the plasmon dispersion at long wave- 
lengths and on the Drude weight, in the same Section 
we briefly discuss interaction corrections to the spin Hall 
conductivity and the renormalization of the spin-orbit 
splitting of the bands due to electron-electron interac- 
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tions. In Sect. IV we present our main numerical results, 
while in Sect.|V|we summarize our findings and draw our 
main conclusions. 



II. GENERAL THEORY 



A. Model Hamiltonian 



where S is the area of the system and 7^ = ^k(0k) is 
given by 



The map 



tan7fc 



a cos (Ok) + f3 sin (Oh) 
P cos (Ok) +asin(6> fc ) 



k -> n cq (fc) = (cos (7/;.), - sin (7*.)) 



(9) 



(10) 



We consider the following model Hamiltonian for a 2D 
electron gas (2DEG), 



soc 



rti. 



(1) 



incorporating the usual parabolic-band kinetic-energy 
term, SOC, and electron-electron interactions. More pre- 
cisely, the first term in Eq. (JTJ) is given by 



#0 =^2 £ ( k )4>l,i 1 PK 

k.i 



(2) 



with i ='|~, I a real-spin label and e(k) = h 2 k 2 / (2irih), % 
being the bare electron band mass. For the SOC term we 
choose a simple linear-in-momentum Rashba-Dresselhaus 
moc 



H 



SOC 



(3) 



Here afj and are Pauli matrices, while a and /3 are the 
Rashba and Dressclhaus SOC constants, respectively. Di- 
agonalization of the sum of the first two terms in Eq. ([!]) 
('Ho+'Hsoc) yields two bands (see, for example, Ref . 142]) . 



e x (k) = e(k) + XkT(Ok) , 



(4) 



with A = ±1 the so-called "chirality" index, Ok the angle 
between k and the x axis, and 



T(0) = yja 2 + /3 2 + 4a(3 sm(0) cos(0) 



(5) 



The Fermi wave vectors for the two bands can be ex- 
pressed in terms of Ok and of the Fermi energy ep : 



, (0) _ , m h T(0k 




m b r( 



1 2 



+ — T2— ■ (6) 



Note that for zero Fermi energy the Fermi contour of the 
minority A = + band shrinks into a single point (i.e. 
0). For any £p > the electron density n can be 



„(o) 



expressed in terms of the Fermi energy as 



m h e F 
nh 2 



2 2 
a 



P 2 



(7) 



The eigenstates of Hq + W-soc corresponding to the 
eigenvalues Q are given by the product of a plane wave 
and a spinor, 



^k,x(r) 



Vs V2\ Ae ~ i7fe 



(8) 



between momentum and the unit vector n cq , which 
parametrizes the noninteracting orientation of the spin 
texture in momentum space, establishes the helical na- 
ture of the model. 

Electron-electron interactions in Eq. |T]) are described 
by the usual two-body spin-independent Hamiltonian 



%int — 



EEE 
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where v q = 27re 2 /(eq) is the 2D Fourier transform of the 
Coulomb interaction (e being a high-frequency dielectric 
constant which depends on the specific semiconductor 
hcterojunction in which the 2DEG is created). This spe- 
cific form of interaction potential applies to a strictly 2D 
system. The finite width of the quantum well hosting 
the 2DEG can be easily taken into account by introduc- 
ing a form factor F(q), which renormalizes the Fourier 
transform v q — > V q — v q F(q) (see, for example, Ref. 43). 

As is common in electron-gas theory^, the electron 
density n will be expressed below in terms of the more 
convenient dimensionless Wigner-Seitz parameter r s : 



1 



irnat 



(12) 



where = eh 2 / (m^e 2 ) is the material Bohr radius. 



B. Equations of motion and plasmons 

Collective modes are isolated poles in appropriate dy- 
namical susceptibilities. Plasmons, in particular, are iso- 
lated poles in the dynamical density-density response 
function x pp (q,aj). Note that we are deliberately not 
denoting the density-density response function by the 
symbol Xpp( t l^ UJ )^ *- e - we are assuming that the sys- 
tem we are interested in is rotationally invariant and 
thus its density-density response function depends only 
on q = \q\. This happens when ft or a is equal to zero. 
The case a = ±/3 deserves special attention and will be 
discussed at a greater length below (sec Sect. IIIF). 

In full generality, this response function can be written 

as 



x P p(q,u) 



1 - V qXpp(<l,v) 



(13) 



where Xppili^) is the s o-called "proper" density-density 
response functiori 20 ! 44 !, which physically describes the 
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density response to the screened potential. The plasmon 
mode can be found by solving the equation, 



and it obeys the standard Heisenberg equation of motion 
(h = 1 from now on) 



1 - v q Xpp(q,u) = 



(14) 



idtPq = \pq,U] = qj 



(p) 



(21) 



In this article we are not interested in the dispersion of 
the plasmon at finite q but only in its limit for q — > 0. 
In this limit we can negleclP^ the distinction between the 
proper and the full causal response function Xpp(9; w )- 
In Sect. |IIC] we prove that 



which is simply the quantum mechanical version of the 
continuity equation. Here the so-called paramagnetic 
current-density operator^ has the following transparent 
form: 



lim lim 3?e Xpp(<7> w ) 

UJ— >0 >0 



V 



(15) 



Jq,x — Yh-q,i 



?x/ 2 



m b 



+ {fi&% - ™\) (22) 



where the quantity T> depends on density and on SOC. 



along the x direction and 



Note the order of limits in Eq. ( 15 ) and everywhere below: 
the limit u> — >• is taken in the "dynamical" sense^, i.e. 
lo 3> wf,a<Z, where wp.A is the Fermi velocity for each 
chiral band. 



5(p) 

■'1,V 



ipk,i + (acj x q - , (23) 



Before proceeding with the formal proof of Eq. ( 15 1 we duced the spin-density operators 



highlight its main physical consequences. Using Eq. ( 15 ) 



order in q, 



in Eq. (14) and solving for ui we find that, to leading 



u p \(q -> 0) = 



*V2 = 



/ 2nne 2 
m pl e 



,V2 



where we have introduced the plasmon mass 



V 



(16) 



(17) 



The plasmon mass is thus completely controlled by the 
quantity T>. In the same limit the imaginary-part of the 
low-frequency a.c. conductivity cr(uj) = ie 2 u)Xpp{u)/q 2 
has the form 



3to a(u>) —} 



V 



(18) 



The a.c. conductivity is a causal response function, which 
implies that its poles can only lie infinitesimally below 
the real-frequency axis, i.e. a(cu — > 0) oc (cv + irj)^ 1 . It 
then follows that the real-part of the conductivity has a 
(5-function Drude peak at u> = 0: 



5Re ct(w) = VS(u) . 



(19) 



The quantity T> introduced in Eq. (151 is thus precisely 
the Drude weight. In the presence of disorder the 5- 



function peak in Eq. ( 19 1 is broadened into a Drude peak, 



but the Drude weight is preserved. 



C. Rigorous definition of the Drude weight 



We now proceed to demonstrate Eq. (151 using the 
"equations-of-motion" approach. The density operator 
corresponding to the Hamiltonian ([T]) is given by the 
usual expression 



Pi 



h.i 



(20) 



along the y direction. In Eqs. (|22|)-(|23|) we have intro- 

(24) 



We now introduce the causal linear-response functions 
Xab{u), which are defined by the Kubo "product"^, 



S 



dt([A{t),B(0)])e lu,t e 



iut — 7]t 



(25) 



where the symbol (O) denotes the expectation value of 
the operator O over the exact interacting ground state 
and r\ — > + is a positive infinitesimal. The dynamical 
response function ((A;B)) U obeys the following identity 



{{A;B)) U = -([A,B] 



-((d t A;B)) u , 



((A;B)) U 



-([A,B])-U(A;d t B)) u 



(26) 



(27) 



Using the continuity equation (plj) and Eqs. (|26 



(27), the density-density response function Xpp(?i w ) can 



be expressed in terms of the longitudinal paramagnetic 
current-current response function as, 



S 
1 1 

S lo 



((q-jq Ph :P-q))^ 



lq-([j q p) ,p- q ]) , K(q-j q P) ;q-j^)), 
S 



;(p)\ 



U3* 



S 



(28) 



We remind the reader that in the presence of a vector 
potential A\~ the physical current-density operator j q is 
related to the paramagnetic one by 

^=^ + m ^£^ fe/ 3 fc . (29) 
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The paramagnetic current-current response function 



1 



(30) 



(here £, to label the coordinate indices) is thus related to 
the physical one by the simple equation^! 

Ti 

Xjdm («?> w ) = hm + XAp) ,-Cp) (g, w) . (31) 

TOb J < Jm 

For generic values (a,/3) of the SOC constants, the 
dynamical response functions of the model described by 
Eq. ([I]) are anisotropic, i. e. they depend on the direction 
of q. However, in the cases of pure Rashba (j3 — 0) or 
pure Dresselhaus (a — 0) SOC the ground state of the 
Hamiltonian Hq + Hsoc is rotationally invariant: for the 
sake of simplicity, in what follows we will restrict our at- 
tention to these two "extreme" cases. From now on we 
assume a ^ and j3 — 0. In Sect. |IIIF| we will com- 
ment on how our results change in the pure Dresselhaus 
(a = and (3^0) case and in the special case a — ±(3. 
Last but not least, we also assume to be in the regime 
in which both chiral bands are occupied (e-p > 0). In 
this situation the Fermi surface consists of two concen- 
tric circles. At low enough densities the topology of the 
Fermi surface changes dramatically, the occupied states 
becoming an annulus in momentum space. We will not 
tackle interaction effects in this interesting but hard to 
achieve experimentally regime. 

In a homogeneous and isotropic liquid we can de- 
compose the tensor Xjdmili 1 ^) its longitudinal and 
transverse components with respect to the direction of q: 

I \ I \9iQm 



Xt(<7,w) 5t 



(32) 



Using this definition we immediately end up with the 
following result 



Xpp{<l,u) = 



l q-([j q P \p- q ]) 
S uj 2 



q 



m b 



(33) 

We stress that Eq. ( 33 ) is exact (provided that the ground 
state is homogenous and isotropic). 

The commutator on the r.h.s. of Eq. p3| can be cal- 
culated easily: indeed, the portion of the paramagnetic 
current operator due to spin-orbit coupling [second terms 
on the r.h.s. of Eqs. ( 22 )- ( 23 )] is proportional to the spin 



operator & q only, which commutes with the density oper- 
ator p^ q . Thus the commutator is found to be equivalent 
to that of the 2DEG without any spin orbit coupling. It 
is related to the so-called f-sum rule^ and is given by 

1 



S 



TO b 



(34) 



Using Eq. ( 34 ) into Eq. ( 33 ) we are left with the following 



crucially important relation: 



X P p{Q,<- 



Xh(q,oj) 



(35) 



Note that the f-sum rule is crucial for the cancellation 
of the diamagnetic n/rrib term in the square brackets on 
the r.h.s of Eq. ([33]). 

Eq. (35) is identical in form with Eq. (15) provided 



that we identify T> with the following dynamical limit: 



V 



lim lim 5Re xl(<7, w) 



w->0 <j->0 



(36) 



This equation is extremely important because it gives us 
an operational definition of the Drude weight. In order 
to calculate it we need to compute the dynamical limit of 
the real part of the longitudinal current-current response 
function xl(<7, w). Such a microscopic calculation will be 
carried out below in Sect. IIIII within the so-called time- 
dependent Hartree-Fock approximation. 



D. Broken Galileian invariance 

In a standard 2DEG without spin-orbit coupling (a = 
(5 = 0) the longitudinal current-current response function 
obeys the exact relation 



lim limXL(?,w) 

Ul— fO q— ^0 



TO b 



(37) 



a nonperturbative result (i.e. valid for any strength of 
electron-electron interactions as long as the 2DEG re- 
mains in a translationally-invariant and homogeneous 
ground state), which is completely independent of com- 
plicated exchange and correlation effects. In this case the 
Drude weight becomes T> = irne 2 /m^ and the plasmon 
mass reduces to the bare electron mass, m p i = rn b . 



The physical reason behind the exact result (37) is the 
following. In the limit q — > xl(9>w) measures the re- 
sponse of the system to a homogeneous time-dependent 
vector potential A(t), i.e. to a homogeneous electric 
field E(t) = ~c~ 1 dA(t)/dt. In a system with a single 
parabolic band the usual replacement p —> p + eA(t)/c 
implies that a uniform vector potential couples identi- 
cally to all the electrons and thus only to the center-of- 
mass motion. This is immediately seen in first quantiza- 
tion: 



H (A) = 



1 

2mb 



Ho 



Pl + -A(t) 



3 cm- A(t) + 0(A 2 ) , (38) 



where Pcm = Yli Pi 1S ^ ne centre-of-mass momentum. In 
the last equality terms of order A 2 have been neglected 
since we are interested in the linear-response regime. 
Electron-electron interactions are thus completely trans- 
parent to A(t), since the latter does not probe the relative 
motion of electrons. 

Eq. ( p7| can be derived by a classical Newton's equa- 
tion for the ccntre-of-mass coordinate Rcm'- 



TO b TV 



d 2 R, 



CM 



dt 2 



-eNE(t) 



e N dA(t) 
c dt 



(39) 



5 



where N is the total number of electrons. Integrating 
this equation we find VcM.(t) = [e/(mbc)] A{t) or 



4 P J (*) = nV C u(t) = — G A{t) 
y mb c 



(40) 



i.e. Eq. (37) 



To see more formally why Eq. (37) comes about, we 
can use the exact-eigenstate (Lehmann) representatiorPS 
for the current-current response function: 

\0\jS\n){n\j^ljO) 



Xjdmil'V) = Sir, 

m b 

;(p) 



S 



10 - UJ n0 + ITj 



m V jJn)(n\j^\Q) 
to + 0J n0 + ir\ 



(41) 



where the limit rj — > + is understood. In a translation- 
ally invariant system, with or without SOC, the exact 
eigenstates \n) are eigenstates of the total momentum. 

In the absence of SOC, moreover, j) „— coincides with 
the total momentum [see Eqs. (22)-(23)] and thus for 
a — ft — and q — > the second term in Eq. (41 ) van- 



ishes and one is left with Eq. (|37|) 

•(p) 



In the presence of 



SOC, however, jfq =Q does not coincide with the total 
momentum and thus Eq. (37) ceases to be true. 



When a (or /3) is non-zero we have 



lim limxL(<?,w) ^ — 
w->0(3->o mb 



(42) 



Deviations from the trivial n/rrih result are due to both 
single- and many-particle effects^. 

The single-particle contribution to the long- wavelength 
low-energy limit of xl(9,w) can be found quite easily. In 
Sect. MI Cl we will show that if electron-electron interac- 
tions are neglected 



lim lim Xl } (?,w) 

uj— q— J-0 



2^0 



m b 



(43) 



where vq — /tt is the usual 2D parabolic-band density- 
of-states in the absence of SOC. The rest of the paper is 
mainly devoted to quantifying interaction-corrections to 
Eq. pi). 



E. Failure of the random phase approximation 

Before concluding this Section, we would like to em- 
phasize that the popular random phase approximation 
(RPA) is not capable of capturing the subtle renormal- 
izations of the Drude weight due to many-body effects. 

By definition, within RPA the proper density-density 
response function Xppil: ^-0 is approximated with its non- 
interacting valued 



~ / ^ RPA (o) 

Xp P {q,u) x P p 



1 ,.(o) 



Xri*,") ■ ( 44 ) 



When Eq. ( 44 ) is substituted in Eq. ( 36 1 one finds im 



mediately that the RPA Drude weight is identical to its 
nonintcracting value: 



£> RPA = lim -limx[° \q,u)) 

uj— *Q q—tQ 



2^0 



m b 



(45) 



More physically, the reason why RPA does not capture 
the subtle interaction renormalizations of T> is the follow- 
ing. During a plasmon oscillation the Fermi circle oscil- 
lates back and forth in momentum space. Due to SOC 
this oscillatory motion of charge excites spin oscillations. 
Exchange interactions are of course very sensitive to the 
spin degrees-of-freedom. The RPA, however, is simply 
a time-dependent Hartree theory^, which treats exactly 
only the self-consistent electrical potential, 



(46) 



created by the electrons displaced away from the equi- 
librium position in the presence of the neutralizing back- 
ground, while completely neglecting the self-consistent 
exchange field associated with the spin degrees-of- 
freedom. From this argument it clearly emerges that the 
minimal theory which can capture interaction-corrections 



to Eq. ( 43 ) is the time-dependent Hartree-Fock theory. 



III. MICROSCOPIC TIME-DEPENDENT 
HARTREE-FOCK THEORY 

In this Section we present a microscopic theory of T> 
that takes into account electron- electron interactions in 
an approximate manner. As discussed in Sect. II E[ the 
minimal approximation that captures the renormaliza- 
tion of V due to many-body effects is the so-called time- 
dependent Hartree-Fock approximation (TDHFA). One 
of the pleasant properties of the TDHFA is that it is 
exact to first order in Coulomb interactions. Other ad- 
vantages, such as its relative simplicity, will be evident 
below. 

As we have amply discussed in the previous Sections, 
we want to study the response of the system described 
by the Hamiltonian |l]) to a weak homogeneous external 
time-dependent electric field directed along, say, x. In 
the gauge in which the scalar potential is zero the elec- 
tric field is simply described by a time-dependent vector 
potential: E{t) = — [c~ 1 dA(t)/dt] x. The vector po- 
tential enters the Hamiltonian ([T ) via the usual minimal 
coupling p — > p + eA(t)x/c. The parabolic-band part 
becomes 



Ho(t) = 

h . i 



-A(t) 
c 



-i>l,iTpk,i , (47) 
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while the SOC part reads 



Hsoc(i) = a {<i k y ~ 4 { k * + - c A (t)\ } $k,j ■ 

(48) 

Neglecting terms C(A 2 ), which are beyond linear- 
response theory, we can write the sum of the two terms 



in Eqs. (47|-(48| as 



Ha(t)+H SO c(t) = H + H SO c + PcuMt) 



(49) 



Thus, due to SOC, a magneto-electric effect appears^ 
a uniform electric field applied along the x direction acts 
as a uniform magnetic field in the y direction [last term 



in the r.h.s. of Eq. C91]. Here df ot = a v 



'q=0- 



Electron-electron interactions are treated within the 
Hartree-Fock (HF) mean-field theory in which the two- 
body term in Eq. (11) is approximated as^ 



(50) 

where (...) (: ... :) denote the expectation value over 
(normal ordering with respect to) the HF ground stated 
At this point we introduce the spin-density matrix, 



Wk^ J k',j) = S k ,k'Pij(k) , 



(51) 



which just assumes that the mean-field ground state is 
translationally invariant. The interaction contribution 
to the total Hamiltonian reads 



Hint = —a^2^2 v k-k'Pji(k') : i>l^ k „ 



(52) 



k,k f i,j 



We parametrize the spin-density matrix Pij(k) in a com- 
pact formal in terms of the occupation factors, rik,±, of 
the noninteracting Hamiltonian Hq + Hsoc m the eigen- 
state representation and in the absence of A(t): 



pij(k) 



n k ,+ + n k - 



rik,+ - n k - 



h{k) ■ a-ji . (53) 



Here h(k) is a unit vector on the 2D plane which de- 
notes the orientation of the spins in the total "effective" 
magnetic field. The idea behind this parametrization is 
that a homogeneous external field (a field with q = 0) 
cannot change anything but the orientation of the spin, 
which is encoded in the unit vector h(k). Note that in 



the absence of the external field n(fc) = n cq (fc). Eq. (53) 



is nevertheless approximate since it assumes the absence 
of interaction effects in the ground state of the system. 
More explicitly, the Fermi wave vectors are renormalized 

. We will 



by electron-electron interactions^, fcp + — » fc p.± 
come back to this point below in Sectr |lHA"} 



Using Eq. (53 1 in Eq. (52 1, the total mean- field HF 



Hamiltonian can be written as 

Hhf = : 4>l,i feSo(fe) + tra ■ B(k)} j, hJ : (54) 

k,i,j 

where the HF fields are defined by 

B (k) = E(k) + —P* M A(t)- [ ~v h -kif + {k') 

(55) 



and 



B(k) = h(k)-J ^ Vk _ k ,f_( k ')h(k') 



In Eq. (56), 



h(k) = —a-A(t)y + ak n oq (fc) 



(56) 



(57) 



is an effective magnetic field, which has the external 
magneto-electric component!^ of modulus 



■Bcxt = a " A 
c 



(58) 



arising from the external vector potential and an internal 
component oc n cq (fe), while the last term is the exchange 
field due to the electron-electron interactions with 



and 



f±(k) 



n(k) 



n k ,+ ± n fc ._ 



B(k) 



\B(k)\ 



(59) 



(60) 



The noninteracting band-eigenstate occupation factors 
are given by 



n k ,± = 6(ef - e±(fe)) 



(61) 



where Q(x) is the standard Heaviside step function. As 
we have already emphasized above, for pure Rasha SOC 
(P = 0) the momentum occupation factors n k ± are 
rotationally-invariant and depend only on k = \k\ (the 
same is true also for pure Dresselhaus SOC, a = 0). It is 
thus extremely convenient to decompose the spherically 
symmetric inter-electron interaction v k - k i in angular mo- 
mentum components, 

Vk-k- = £ V m (k,k')e im ^- e ^ , (62) 



771 — — OO 



with 



V m (k,k') = 



2tt 



d6 



V l\q=\k-k'\ 



(63) 



9 being the angle between k and k! . 
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A. Equilibrium HF theory 

In the absence of the external electric field, i.e. A(t) = 
0, the unit vector n(fc) coincides with the equilibrium 
one: 

h{k) n cq (fe) = (sin (9 k ), - cos (0 fc )) . (64) 



Substituting Eq. ( 62 ) in Eq. ( 56 ) and performing the an- 



gular integration over 9k> , we find that the equilibrium 
solution of Eq. ( 56 1 reads 



B cq {k) 



xk 



dk' 

2n~ 



k'f-(k')Vi(k,k') 



n cq (k) 



(65) 

As expected, in the absence of the electric field, B(k) is 
oriented along n cq (k) and it is isotropic. The modulus 
of B cq (k) is simply 



|B oq (fe)| =afc + Ex(fc) , 



(66) 



and depends only on k with the self-energy £i(fc) defined 
by 



dk' 

~2n~ 



k! /_(fc')Vi(fc,fc') 



(67) 



For Ep > the factor /_ in the integrand, being the 
difference in the occupation of the two bands, picks up 
contributions only from wave vectors fc'_in the interval 



[k F ' + , k F _], where k F ' x is given by Eq. (|6|) with j3 = 0. 
The self energy can thus be written as 



Ei(*) 



1 

47T J k («) 



dk' k'Vi(k,k') 



(68) 



In a completely analogous manner, it is possible to find 
the equilibrium solution of Eq. (551 which reads 



B ,eq(k) = e{k) 



dk' 
2tt 



k'f + (k')V (k,k') , (69) 



or, B ,eq(fc) = e(fc) + So(fe) with 



1 

~2^ 
1 

47T L(°) 



dk' k'V (k,k') 



dk' k'V (k,k') 



F,+ 



Finally, the complete HF bands are given by 
£W,A(fc) = s(k) + S (fc) + X[ak + S^fc)] , 



(70) 



(71) 



and the quasiparticle effective mass m\ for the A-th band 
can be defined as 



4 0) a _ dE HFA (k) 



Ok 



(72) 



Rigorously speaking, the HF bands and the corre- 
sponding Fermi wave vectors (&f,±) should be calculated 
in a fully self-consistent manner. We have done this and 
we find that the repopulation of the energy bands in the 
ground state due to interactions is a very small effect. In 
fact, the difference between fcp ^ and kp t ± is less than 
0.5% over the entire range of parameters we have consid- 
ered. We have thus ignored this small effect throughout 
this article and used fc F °' ± in all calculations. 



B. The non-equilibrium problem: linearization of 
the HF equation 



We now proceed to solve Eq. ( 56 1 in the presence of the 



external electric field by linearizing it around the equi- 



librium solution B cq (k). To this end we write 



B(k) = B cq (k) + SB{k) 



(73) 



and 



h{k) = n cq (k) + —^-l + 0{{8Bf) , (74) 

where SB±(k) — SB(k) — n cq (k)[n cq (k) ■ 5B(k)} is the 
component of SB(k) perpendicular to h cq (k). We now 
make the following Ansatz for SB(k): 

SB{k) = [5B htl {k) cos (0 fe )]n eq (fc) 

[SB TA (k)sin{9 k )}zxn cq (k) . (75) 

We note that the Ansatz has to be consistent with the un- 
derlying model Hamiltonian. Indeed, for the pure Rashba 
model, n cq (k) = k x z and z x n cq (k) = k. Since the 
magneto-electric field is in the y-direction [see Eq. ( 57 )] 
its components along z x n cq and n oq are proportional to 
sin(6*fc) and cos(#fc), respectively. This justifies the par- 
ticular form of Eq. ( |75| . Using Eq. (75) in Eq. ([74]) we 
find that 

8Mk) = m-n^) = - S ^g^ k. (76) 



Substituting Eqs. (75) and (76) in Eqs. (56), integrat 



ing over 9k', and keeping only terms that are linear in 
5B(k), we find 



\[5B ul {k) 



5B TA {k) 
dk' 

47T 



— -E>cxt 

k'f.(k')v (k,k') SB ^ ik,) 



■Bcq(fc') 



(77) 



and 
1 



6B TA (k)} = 

SB TA {k') 
\B cq (k')\ ■ 



dk' 
An 



k'f.(k')V 2 {k,k') 



(78) 
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Summing and subtracting these two equations we fi- 
nally find the following integral equation for the trans- 
verse <5-Bt,i(&) component: 



/>oo 

6B TA (k) = B cxt - / dk' K. T {k,k')SB TA (k' 
Jo 

where the kernel K.T{k, k') is defined by 
JC T (k,k') 



(79) 



We need to evaluate the longitudinal response xl to 
a uniform vector potential A in the to — > limit. We 
thus have to evaluate the change in longitudinal physical 
current due to a uniform electric field applied, say, along 
the x direction: 



Sj x = lim 

A->0 



(jq=0,x)A — (jq=O,x)A=0 



(87) 



1 V (k,k') + V 2 (k,k>) 

4n IAk) \B cq (k>)\ 



1 k'f _ ( ,j ) Mk,k') + V 2 (k,k') 
4tt 7 1 1 ak' + Etik') 



• (80) 



where the physical current operator j q has been intro- 
duced in Eq. (29). Recalling the definition of the spin- 
density matrix Pij(k) we find that 



Once Eq. ([79]) has been solved self-consistently for 5 3x = ^ ^y + ^ ~^ 6pit ^ ~ a ~S a ij 5p v ^ 

h.i 



SB^^i(k)^ the longitudinal component SB^_i(k) can be 
calculated from 



5B LA (k) = B cxt - / dk' JC L (k,k') 6B T .i(k') 
Jo 



with 



KAk k')-U'f 

^L(fc,fc)- 47r fc/-(fc) afc , + Si(fc0 



(81) 



(82) 



For future reference, it is very convenient to rewrite 
Eq. (79) in a dimensionless form. To this end, we scale 



all the wave vectors with the 2DEG Fermi wave vec- 
tor kp = \[2wn in the absence of SOC, all energies 
with ep.o = fcp/(2rrib), the pseudopotentials V m with 
2ire 2 /(efcp)j and, finally, we introduce the dimension- 
less SOC Constantsa = rriba/k-p and the dimensionless 
quantity u = SB^.i/ B cxt - From now on, symbols with a 
bar over them denote dimensionless quantities. In these 
units Eq. ( 79 ) reads 



( , , , r 3 f A - , , V (x,x') + V2(x,x') , 

u[x) = 1 H = / ax x ^— — u(x ) 

V ; 2V2 J A+ 2ax' + H 1 (x') v ; 



Here x = k/kp, x' = k' /kp, 



A- 



k 



(<>) 

F,± 



ky, 



(83) 



(84) 



and Si(x) is the dimensionless version of the self-energy 
introduced in Eq. (|68|: 



£i(a;) 



£f,o 



^F(x) 



with 



F(x) 



dx' x'Vi(x, x') . 



(85) 



(86) 



C. Interaction corrections to the Drude weight and 
renormalization of the in-plane spin susceptibility 

We are now in the position to evaluate the interaction 
corrections to the Drude weight from the definition in 
Eq. ((361). 



k,i,j 



Since the diagonal components of the spin-density ma- 
trix do not change under the application of a uniform 
magnetic field we get the following important relation 



Sjx 



A 



mh c 



(89) 



where 
8<j v 



I^4^ i( fe) = 2/^/_(fc) 6n y 



= -2 



(27 

d 2 k . a 

r/-(fc) ,- — sm 



9 k ) ■ (90) 



(27r) 2 " v ' \B eq (k) 
Performing the angular integration we finally find that 

S* v = X*»avBe*, (91) 

where the magneto-electric field B cxt has been introduced 
above in Eq. ( |58[ ) and where the in-plane spin suscepti- 
bility is given by 



1 



X<y y <y y 



B c: 

2 -/A 



dk 6B TA (k) 
4^ "|B cq (fe)| 



A_ 



dx 



xu{x) 



2ax + Si(x) 



(92) 



In the noninteracting r s —¥ limit the vertex correction 
u tends to unity and the self-energy to zero: in this 
limit Eq. (92) reproduces the well known result for the 



in-plane spin susceptibility of a 2DEG with Rashba SOC, 

(0) /o 

i.e. xJ<jy = v ol 2 - 



Finally, using Eq. ( 36 ) we find that the Drude weight 
is given by 



V = Tie/ 



Sjx 

eA/c 



a x<yy<jy 



(93) 



This is the most important result of this work. It states 
that the corrections due to SOC and many-body effects 
to the universal Tre 2 n/m^ Drude weight of a standard 
parabolic-band 2DEG are completely controlled by the 
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uniform in-plane spin susceptibility x<y y <y y in the dynam- 
ical limit. Note that, even though the particular deriva- 
tion we have given in this Section seems to be related to 
(and thus dependent on) the TDHFA, Eq. (93 1 is exact 
and stems directly from Eq. ( 22 ) . 

In the noninteracting limit Xa y a y ^o/2 and thus, in 
the same limit, 



Vn 



ire 



n 



2^0 

a y 



(94) 



In Sect. IV we will present numerical results for the ratio 
T>/T>q as evaluated from the HF expression for the in- 
plane spin susceptibility in Eq. ( 92 ) . Normally electron- 
electron interactions enhance the spin susceptibility: we 
thus anticipate that the Drude weight of the interacting 
system is smaller than its value T>q in the absence of 
interactions. 

Before concluding this Section we derive a semi- 
analytical expression for x<y y a y U P to first order in the 
coupling constant e 2 . To this order of perturbation the- 
ory the solution of Eq. (83) can be found analytically 
with the result 



u(x) = 1 + 



(95) 



where 



„(*)= [ A - dx >M^±M^l . (96) 

JA + 2 



We notice that the perturbative solution ( 95 ) is not of the 
first order in r s , since a = m\,oi/k-F and A± also depend 
on density via the Fermi wave vector. In the presence 
of SOC, interaction effects are not solely controlled by 
r s . Substituting Eq. (951 in Eq. ( 92 ) and expanding the 



ratio in the integrand of this equation in powers of e up 
to first order we finally find that 



Xa y a y 



= 1 



r s A 



with 



A = 



dx 



F(x) 



(97) 



(98) 



A plot of A as a function of r s for different values of 
a is reported in Fig. [I] we clearly see that A is positive 
and that thus the in-plane spin susceptibility is enhanced 
by electron-electron interactions (at least to first order in 
e 2 ). In the high-density and/or weak-SOC limit (nibQ "C 
fep or, equivalently, a < 1) we can approximate A in the 
following manner: 



0.30 
0.25 
0.20 
^ 0.15 
0.10 
0.05 



0.00 



a = 5 x 10 11 eV m 
a = 10" 11 eV m 



0.0 0.2 0.4 0.6 0.. 



1.0 1.2 1.4 



FIG. 1: (Color online) The quantity A in Eq. (981 as a func- 
tion of r„ for two values of the Rashba SOC strength a. Notice 
that A > and thus x<? y <y y IXalav > !• 



In the same limit 



g(l)-F(l) -> 2a 



F (l,l)+Va(l,l) 



4a 
37 ' 



(100) 



the last equality being valid only for unscreened Coulomb 
interactions [see Eq. (1211 below]. In this case and for 
ol — ~> we find A — > 8a 2 /(37r). Using this result in 



Eq. (97) we find a rigorous result for the spin suscep- 



tibility enhancement to linear order in r s : 



Xa y a y 

A(jV(jV 



V2 



(101) 



We finally remark that, in the oversimplified case of 
ultra-short-range interactions, 



constant 



27re 2 



(102) 



Eqs. ( 83 ) and ( 92 ) can be solved analytically. In this case 



indeed all the moments V m (k, k') of the inter- particle in- 
teraction but the m = one are zero. The solution of the 
integral equation (83) is a constant u = [1 — ^kob) -1 ] -1 



and the in-plane spin susceptibility turns out to be equal 
to u: 



XaVa-y 



1 



1 - (2kob) 



> 1 



(103) 



D. Interaction corrections to the optical spin Hall 
conductivity 



A 



(A--A+) 



F(x) 



2%(1)-F(1)] 



(99) 



It turns out that the in-plane spin susceptibility x& y <? y 
introduced in the previous Section controls also the "opti- 
cal spin Hall conductivity" asn(u>) of the Rashba model. 
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This was first shown by DimitrovcP^I. For the sake of 
completeness, we briefly summarize here the key steps of 
the derivation. 

The spin operator at q = 0, a^ ot , satisfies a simple 
equation of motion: 

id t *l t = Kt,n] = [<J v tot ,Hsoc] = -4im b af y , (104) 

where the q = f-spin current operator in the y direc- 
tion, jy, is defined by (see for example Ref. 

K 
m b 



2 ^ 



(105) 



The spin Hall conductivity ctsh(w) describes a z- 
polarized spin current flowing in the y direction in re- 
sponse to a homogeneous (q — 0) electric field E = E x x 
along the x direction: 



(106) 



From Eq. ( 104) it is immediately evident that in the d.c. 
limit ujt — > (r is the electron-impurity scattering time) 
the spin Hall conductivity is zero since in a steady state 
(<9t<5f ot ) = 0. This is the limit that is relevant to d.c. 
transport. The vanishing of the transport spin Hall con- 
ductivity in the Rashba model has been widely discussed 
in the literature (see e.g. Ref. |4"9"|) . 

Here we are interested in the high-frequency or clean 
limit, ujt — > oo, which can in principle be probed in time- 
resolved experiments with photo-excited carriers^ or in 
ballistic transport. In this limit the following analysis is 
particularly useful. We first notice from Eq. ( 106 1 that 



the optical spin Hall conductivity is related to the spin- 
current response function by 



<tsh(w) 



where ji P ^ — j^o x 1S the x component of the q = 



X 3 ,, P ,H = -((j-j^)) w , (107) 



;(P) 



paramagnetic current operator in Eq. (22). We then sub- 



stitute Eq. (104) in the response function on the r.h.s. of 



Eq. (107) and we use Eq. (26). We get 



«^;j! p) » 



1 



4mbO 
^< 

, U°totJ Jx 



((3 t of ot ;i' 



(p) 



4m b 



tot! "tot//" 1 



(108) 



where we have used that ((it^;?^))^ = and that 
[<5f ot , Pq M ] = 0. The former is a consequence of the fact 
that total momentum is a conserved quantity (in the ab- 
sence of impurities) even in the presence of Rashba SOC, 
while the latter is a trivial commutation rule. Using 
Eq. |l08l ) in Eq. ( fT07| we finally find 



csh(w) 



4mi 



(109) 



In the high-frequency or clean lot — > oo limit and 
for noninteracting electrons we have x<y y <yy ~^ v^jl 
and thus Eq. (109) gives the well-known "universal" 
value <7sh(wt — > oo) = e/(87r). As we have seen 
above, however, electron-electron interactions enhance 
the high-frequency spin susceptibility, thus yielding an 
enhancement of the optical spin Hall conductivity. Us- 
ing Eq. ( 101 ) we immediately find that for lot — > oo 



csh 




(110) 



Eq. (110) has to be compared with Eq. (36) in Ref. [27] 



where an identical result was found modulo the sign of 
the second term in round brackets. Dimitrova indeed 
predicts a suppression^ of the spin Hall conductivity due 
to interactions rather than an enhancement. 

For the case of ultra-short-range interactions [see 



Eq. ( 102 )] we find that 



- — 11 — 



> 



8?r 



(111) 



Before concluding this Section we would like to men- 
tion th at it is possible to derive a relation similar to that 
m Eq. ( flQ9| for the spin Galvanic effect^, i.e. the gener- 
ation of a charge current in the x direction in response to 
a homogeneous Zeeman magnetic field B = B y y applied 
along the y direction: 



jx = <y SG B y . 



(112) 



Following a procedure analogous to the one that led to 
Eq. (fl09|), we find 



(113) 



where g is the material Lande gyromagnetic factor and 
/iB is the Bohr magneton. 



E. Interaction-induced enhancement of the Rashba 

SOC 



From the functional form ( 71 1 of the HF bands of the 



Rashba model it is evident that, for a given density, the 
real part of the a.c. conductivity 



ifte o~(lo) = — e 2 lim 



q-yQ q z 



Xppil, 



(114) 



is finite (i.e. absorption occurs) only in a finite interval 
of frequencies: A + < lo < A_, where 



A+ = EnF,+ (kf;+) - £hf,-(4,+) 



A_ = £W,+ (A^J-£W,-(*# V -) 



(115) 
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In the noninteracting limit these bounds ar d 23 l 24 h = 
(°) or,^ a(°) — 9™j-( ) j n the interacting case we 



2ak%> + and = 2afc F 



can define A_i 



25^fel 0) , and A_ 



2a_fcp°^_, where 



5+ 
a 


= H 


r s A_ 
2y/2a 


r 

Ja+ 


5_ 


= H 


r s A+ 


r 


a 


2y/2a 


7A+ 



das' cc'yi(A + ,a; / ) 



dx' x'Vi(A-,x') 



(116) 



We can thus view a± as effective Rashba SO C stre ngths 
renormalized by electron-electron interactions^^!] xhe 
dependencies of a± on r s and a will be illustrated below 
in Sect. 13 

In high-density and/or weak SOC limit (a -C 1) we 
find that A + = A_ = 2akp with (restoring physical 
dimensions for a moment to make contact with earlier 
work) 



o 

Q 



1 



m b 



dJS 



COS (0) 



,(H7) 



2nh' Jo 2tt ^ig=2fe FSi „(e/2) 
in perfect agreement with the work by Chen and RaiklPS. 

F. The pure Dresselhaus case and the a = ±/3 case 

Before turning to the numerical results, we would like 
to mention that all our results apply equally well to the 
pure Dresselhaus model, i.e. for a = and finite /3. 
Indeed, replacing the Rashba interaction by the Dressel- 
haus interaction has the only effect of changing the phase 
7^ of the cigenspinors ^k,x(r) in Eq. ^ by tt/2, leav- 
ing the band energies in Eq. Q unchanged. We have 
also checked that the statement at the beginning of this 
Section is true by applying the TDHFA to the pure Dres- 
selhaus model. 

The model with a = ±/3 is much more subtle. For 
a = 0, for example, Hsoc reads 

"Hsoc = aJ2^kA k ^ + k y)( a f3 ■ ~ VijWkj 

= teJ2fh,ir7jk+4hj , (us) 

k,i,j 

where (Ty = (cr?- - a\- ) / \/2 a nd fc + = (k x + k y )/\/2. From 
the second line in Eq. ( |118| we immediately see that the 
magneto-electric field generated by the application of a 
uniform vector potential A(t) is 



a-(A x + A y )(l,-l) , 



(119) 



i.e. it is parallel to n cq \ a= „ = (1, — 1)/V2, and does not 
affect the spin orientation. Thus, in the case a = ±/3 
a uniform vector potential does not reorient spins. The 
plasmon mass and the Drude weight are thus completely 
unrenormalized by electron-electron interactions. Of 
course the plasmon dispersion at finite q will be sensi- 
tive to interactions. 



IV. NUMERICAL RESULTS 

We now turn to a presentation of our main numerical 
results. As far as the material parameters are concerned, 
in this article we present results for a 2DEG hosted in 
a InAs quantum well. In this material the bare elec- 
tron mass is mb ~ 0.023 m c , where m c is the electron 
mass in vacuum, and the high-frequency dielectric con- 
stant is e rs 15. The material Bohr radius turns out to be 
ob ~ 348 A. As a consequence, a Wigner-Seitz density 
parameter r s — 1 corresponds to a rather low electron 
density, n ~ 2.6 x 10 10 cm -2 . The SOC strength in InAs 
varies in the ranged a ss (1 — 6) x 10 -11 eV m. 

For the numerical calculations we have used a model 
interaction potential of the form 



2Tre 2 



e(q + ^qT F ) 



(120) 



where <7tf = 2/<zb is the 2D Thomas- Fermi screening 
wave vector in the absence of SOC and £ e [0,1] is a di- 
mensionless control parameter; £ = implies unscreened 
Coulomb interactions while £ = 1 implies Thomas-Fermi 
screened Coulomb interactions. The TDHFA is well 
known to overestimate many-body effects when the un- 
screened Coulomb potential is used. (As we have already 
mentioned earlier, when the unscreened Coulomb poten- 
tial is used the TDHFA is exact to first order in e 2 .) 
On the other hand, when statically screened Thomas- 
Fermi interactions are used many-body effects are typi- 
cally largely underestimated. Thus the spirit of the con- 
trol parameter £ is to provide us with upper and lower 
bounds for the strength of interaction corrections to the 
various observables we present in this Section. 

For £ = the coefficients V m (k,k') of the angular- 
momentum expansion in Eq. ( 62 ) can be calculated ana- 



lytically: in dimcnsionless units these are given by 



V m (x,x') 



dt J m (tx)J m (tx') 
x' m r(m + l/2) 

~ r(m + i)r(i/2) 

x 2 Fi{m+ 1/2, 1/2, m+ l,x' 2 /x 2 ) , 

(121) 

for x > x' . Here J m {z), T(z), and 2Fi(a, b, c, z) are the 
Bessel function, the Euler Gamma function, and the hy- 
pergeometric function, respectively. For x < x' one needs 
to interchange i f> i' in Eq. ( 121 ). For £ ^ the pseu- 



dopotentials V m (x, x') have to be calculated numerically. 

Fig. [2] shows the HF bands £ H f,± (k) in Eq. ^ for 
unscreened Coulomb interactions, while Fig. [3] illustrates 
the HF self-energies So(fc) and £i(fc), defined in Eqs. ( 70 ) 
and (68 1, respectively. Note that ^o(k) is negative while 
Si(fcjis positive. 

In Figs. |4^l) and b) we present the minority m* + 
and majority m* effective masses as functions of r s for 
Thomas-Fermi screened interactions, as calculated from 
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tions. This result stems from exchange interactions and 



FIG. 2: (Color online) The renormalized Hartree-Fock energy 
bands EuF.±(k) (in units of £f,o) as functions of k (in units 
of fcp) for r s = 1 and a = 5 x 10 -11 eV m for the case 
of unscreened (£ = 0) Coulomb interactions. The dashed 
(blue) line denotes the minority band [E-HF,+ (k)] while the 
solid (red) line denotes the majority band [EuF,-(k)]. The 
thin black lines are the energy bands for the noninteracting 
case. The vertical lines denote the location of ±k^\. 
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FIG. 3: (Color online) The Hartree-Fock self energies So (A:) 
(solid line) and £i(fc) (dashed line) (in units of £f,o) as 
functions of momentum k (in units of fep) for r s — 1 and 
a = 5 x 10~ 1 eV m. The vertical lines denote the location 
of ±fc 



(0) 



Eq. (72). (It is very well know that, to avoid artifacts 
of the HF theory, it is necessary to screen the Coulomb 
interaction to get meaningful results for the quasiparti- 
cle effective mass. The derivative of the HF quasiparticle 
energy indeed diverges at the Fermi surface and hence 
the quasiparticle effective mass vanishes.) In the case 
of no SOC (a — 0) one finds that the quasiparticle ef- 
fective mass is suppressed by electron-electron interac- 



is the dominant effect at weak couplin, 



,43 



i.e. 



for r s !$ 1. 



The suppression of the quasiparticle effective mass for 
a = shown in Figs. |4^,)-|4Jd) extends to larger values of 
r s because of the strong value of the screening parameter 
(£ = 1) we have used. We notice that the impact of SOC 
is opposite in different bands [this is ultimately du e to 

the dependence of the factor in the l.h.s. of Eq. ( 72 ) 
on a] : as we can see in Fig. [4^,) , SOC further suppresses 
the quasiparticle effective mass in the minority band. On 
the other hand, as shown in Fig. |4|d), SOC enhances the 
quasiparticle effective mass in the majority band. 

In the inset to Fig. [4Jd) we plot the quantities Av\ = 
v x ~ ^aIq^o as functions °f & and f° r a fixed value of 
r s (? s — 0.25), where the quasiparticle velocities v\ are 
defined by 



dE 



HF,A 



(k) 



Ok 



(122) 



k=k { ' 



Differently from the effective mass results, we see that 
Ai?£ is practically the same for both A = ± bands and 
that corrections linear in a are absent, in agreement with 
Refs. 1231551 

In Fig. [5] we report the solution <5I?T,i(fc) of the inte- 
gral equation ( 79 ) for both unscreened and screened in- 



teractions. It is important to note that <5-BT,i(fc) in units 
of the bare effective magnetic field B ext = eAa/{hc) is 
larger than unity. Kinks are seen in SB^ i i(k) at k = k^\, 
which are especially visible at £ = 0. We also clearly see 
how the amplitude of SBt i(k) decreases with increasing 

A plot of the in-plane spin susceptibility Xovov in units 
of the noninteracting value xlyy^y = V Q /2 is presented in 
Fig. [6} As expected, the in-plane spin susceptibility is en- 
hanced by electron-electron interactions. Overscreening 



their strength by setting £ = 1 in Eq. (120) substantially 
reduces the enhancements. Note also that increasing 
the SOC strength a the ratio X<j y <? y I X^y <rv increases. 

Fig. [7] shows the most important result of this work, 
i.e. the renormalization of the Drude weight T> due to 
interactions. There we indeed plot the ratio between T> 
and its noninteracting value T>q. Since the spin suscepti- 
bility is enhanced by interactions, T> is suppressed. The 
suppression is quite large within truly first-order pertur- 
bation theory (£ = 0) and increases with increasing a. 
The plasmon mass is thus enhanced by the interactions 
and thus the plasmon frequency is reduced by the com- 
bined effect of SOC and interactions with respect to the 
standard frequency of plasmons in the absence of SOC. 
The enhancement of SOC due to interactions is il- 
whcre we have plotted a±/a as 
116). These two figures refer to 



lustrated in Figs. [8p 
calculated from Eq 
two different values of £. For the sake of comparison, in 
Fig. [9] we have also plot ted t he weak-SOC result by Chen 
and RaikfPS [see Eq. ( |117[ )]. From Fig. [8] we see that 
the enhancement of SOC is pretty large for unscreened 
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FIG. 4: (Color online) The minority [panel a)] and majority [panel b)] quasiparticle effective masses m± (in units of the 
bare mass mb) as functions of r s for different values of a. These results have been obtained by using the definition (721 
and fully-screened Thomas- Fermi interactions (£ = 1). The solid line represents the classic result by JanalP^ our results for 
a = (filled circles) are in excellent agreement with the analytical expression (16) in Ref. 1541 The inset to panel b) illustrates 
Av± = v± — V±\ a _ (in units of vf = kp/m^) as functions of a and for r s = 0.25 (dashed and dotted lines). The solid line 
is the weak-SOC analytical result (74) in Ref. 1251 Notice that our numerical results extend up to a large value of the SOC 
constant since a = 0.7 corresponds to a ~ 38 x 10 -11 eV m. 
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that for a 



5 x 10~ n eV m, k F °^ + /k F 



0.6 and 4°-/ fc F 



1.3. For a = 10~ n eV m, k^ + /k F 
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Coulomb interactions and that it decreases for increasing 
SOC strength. 



A. Taking into account the density dependence of 
the Rashba SOC 

Until now we have treated the Wigner-Seitz parame- 
ter r s (or density) and the Rashba SOC constant a as 
two independent parameters. This is similar in spirit to 



what has been done for decades in the context of tunnel- 
coupled double quantum wells where the single-particle 
symmetric-to-antisymmetric gap AgAs and density have 
been treated as independent parameters (see e.g. Ref. 1551 
and references therein to earlier work). 

In reality, when a single gate voltage is applied to the 
2DEG to change its density (and thus the r s value) the 
asymmetry of the quantum well which hosts the 2DEG 
changes tooSSl. This in turn changes a. In a simple single- 
band model with infinite barriers the SOC strength a is 



14 



1.35 
1.30 
1.25 
! 1.20 

b 



1.15 

1.10 

1.05 

1.00 
0.0 



1 1 1 1 

■ Eq. (101) 

a = 5 x 10 -11 eV m 

- a = 10" 11 eV m 

e = o 

(a) X* 



J I I L 



J L 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



1.35 
1.30 
1.25 
I 1.20 

b 

^ 1.15 
1.10 
1.05 
1.00 



"i 1 1 1 1 1 r 

a = 5 x 10~ n eV m 
- a = 10- 11 eVm 



^ = 1 
(b) 



J I I I I I L 



0. 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIG. 6: (Color online) The in-plane spin susceptibility Xo-bo-h (in units of the noninteracting value, X^va-y = u o/2) as a function 
of r s and for different values of a. Panel a) Results for unscreened Coulomb interactions (£ = 0). Panel b) Results for 
fully-screened Thomas- Fermi interactions (£ = 1). 



1.00 



0.95 



0.90 



0.85 



n i r 

■ IHILJ.NIII1.., 



i i i r 



(a) 

e = o 

RPA 

a = 3 x 10" 11 eV m 
a = 5 x 10" 11 eV m 



0.80 

0.0 0.2 0.4 0.6 



J I I L 



1.0 1.2 1.4 



1.00 



0.95 



0.90 



0.85 



0.80 



1 1 1 1 1 1 1 






(b) 


N 


£ = i 




RPA 




a = 3xl0- n 


eV m 


a = 5 x 10" 11 

i i i 


cV m 
i i i i 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIG. 7: (Color online) The Drude weight D [in units of the noninteracting value, T>q = 7re 2 (n/mb — a 2 fo/2)], calculated from 
Eq. (931, is plotted as a function of r 3 and for different values of a. Panel a) Results for unscreened Coulomb interactions 
(£ = 0). Panel b) Results for fully-screened Thomas- Fermi interactions (£ = 1). 



given bjismzna 



>(E) 



ea so n 



(123) 



with a so = 117 A 2 for bulk InAs. Here we have used that 
the electric field in the well is given by (E) = n^/e where 
the density of the donors has been approximated by 
the density of electrons n. 

In Fig. [lO] we present numerical results for the ratio 
a±/a calculated by taking into account the density de- 
pendence of a according to Eq. ( 123 1. From this plot we 



clearly see that the difference between the effective SOC 
constants 5+ and a_ becomes negligibly small and that 
the ratio a±/a changes by roughly sixty percent when 
density is changed over three orders of magnitude (for 
unscreened Coulomb interactions). The enhancement of 



SOC due to interactions increases with decreasing den- 
sity. Overscreening Coulomb interactions washes out this 
effect yielding a tiny renormalization over the same den- 
sity range. 



Finally, in Fig. fTTJ we show the spin susceptibility en- 



hancement Xo-» a-y / 'Xo-j Uv as a function of r s calculated 
by taking into account the dependence of the bare a on 
density via Eq. ( 123 1. Note that for unscreened Coulomb 
interactions the enhancement can be as large as 20 — 30% 
for n w 10 10 cm~ 2 (recall that in InAs r s — 1 corresponds 
to an electron density w 2.6 x 10 10 cm -2 ). 
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FIG. 8: (Color online) Renormalized SOC coupling strengths a± (in units of the bare value a) for the minority (+) and majority 
(— ) bands as functions of r s and for different values of a. These results refer to unscreened Coulomb interactions (£ = 0). 
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FIG. 9: (Color online) Same as in Fig. [8] but for f = 1. The thin solid (black) line in both panels labels the result of Ref. 1211 
[see Eq. |117l], which is asymptotically exact in the weak SOC limit. 



V. CONCLUSIONS 

In summary, we have studied the long-wavelength 
plasmon dispersion and the Drude weight of a two- 
dimensional electron gas with Rashba spin-orbit cou- 
pling. We have shown that these measurable quantities 
are sensitive to electron-electron interactions due to bro- 
ken Galileian invariance and we have discussed in de- 
tail why the random phase approximation is not capa- 
ble of describing the collective dynamics of these systems 
even at very long wavelengths. We have then presented 
approximate microscopic calculations of these quantities 
based on the so-called time-dependent Hartree-Fock ap- 
proximation. We have found that interactions enhance 
the plasmon mass and suppress the Drude weight. 

These findings can in principle be tested experimen- 
tally by inelastic light scattering, electron energy loss, 



and far-infrared optical-absorption measurements. In- 
elastic light scattering^ nas already been extensively 
used to measu re the plasmon dispersion in GaAs quan- 
tum well j 60 * 6 H Notable deviations from the predictions 
of the random phase approximation have been observed^ 
at finite momentum transfer q and at low densities. We 
hope that similar studies can be performed systemati- 
cally in asymmetric n-doped quantum wells with tunable 
spin-orbit coupling. 

Last but not least, we have also computed quantita- 
tively the renormalization of the Rashba spin-orbit cou- 
pling constant due to electron-electron interactions and 
the interaction corrections to the clean-limit spin Hall 
conductivity. 

In the future we plan to extend these studies to the 
complete spin-orbit-coupling model Hamiltonian ([3| with 
both a and /3 finite. As already stressed in the main body 
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FIG. 10: (Color online) The renormalized SOC constants a± [in units of the noninteracting value a] are plotted as functions of 
the logarithm of density (expressed in units of cm -2 ). In this figure we have taken into account the density dependence of the 
bare SOC coupling a in the simple approximation given in Eq. (123 1. Panel a) Results for unscreened Coulomb interactions 



(£ = 0). Panel b) Results for fully-screened Thomas- Fermi interactions (£ = 1) 
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FIG. 11: (Color online) The spin susceptibility enhancement 

X<r y c y / X^y o-y as a function of the logarithm of density (ex- 
pressed in units of cm -2 ). In this figure we have taken into 
account the density dependence of the ba re SO C coupling a 
in the simple approximation given in Eq. (1231. 



of this article, this complicates things quite a bit since 
the resulting ground state does not possess rotational 
invariance. It is worth exploring also other spin-orbit- 
coupled two-dimensional quantum liquids such as hole 
gases in quantum wells which have a very rich single- 
particle band structure^. 
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